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Abstract 

We study the entanglement emergence in a dipolar-coupled nuclear spin- 1/2 system cooled using 
the adiabatic demagnetization technique. The unexpected behavior of entanglement for the next- 
and next-next-neighbor spins is revealed: entangled states of a spin system appear in two distinct 
temperature and magnetic field regions separated by a zero-concurrence area. The magnetic field 
dependence of the concurrence can have two maximums which positions are determined by the 
initial conditions and number of spins in a chain. 



1 



I. INTRODUCTION 



Prospects for practical application of quantum technology and devices based on entangled 
state stimulate intensive qualitative and quantitative research of quantum entanglement in 
various physical systems. Among these systems, a system of nuclear spins presents excel- 
lent theoretical and experimental models for studying the entanglement properties. Systems 
with a large number of spins with various models of the spin interaction have been described 
in a number of papers, see e.g. jl-5|, and the references therein. However, a major part 
of the performed studies considers the temperature and field dependences without taking 
into consideration the realistic cooling technique for achieving low temperature at which 
entanglement appears. Adiabatic demagnetization (AD) performed by variations of an ex- 
ternal magnetic field is a very useful and effective technique for the attainment of proper low 
temperatures, of the order of microkelvins, at which entanglement appears in spin systems 

d-liil- 

Recently, using adiabatic demagnetization in the rotating frame_(ADRF) 12Ml4J| , en- 



151 ] . In contrast to the 



tanglement in a dipolar-coupled nuclear spin system was studied 
conventional AD method, ADRF uses an external magnetic field which much stronger than 
internal dipolar field and a radiofrequency field with an offset from the Larmor frequency. 

It was shown that there is no entanglement between remote spins, although the considered 
model included dipole interactions between these spins which is the typical condition to 
create entangled states. 

In this paper we consider the entanglement emergence in a one- dimensional dipolar- 
coupled nuclear spin- 1/2 system cooled using the conventional AD technique. We consider 
the case where the Zeeman energy is of the order of or even less than the dipolar interaction 
one. 

The structure of the paper is as follows: in the next section, we describe the Hamiltonian 
for a spin system in an external field, conditions of the AD realization and analyze variations 
of the spin temperature T and magnetization M at AD. Then we consider the pairwise en- 
tanglement. Discussion of the results is given in the final section. The numerical calculation 
of the spin temperature, magnetization, heat capacity and concurrence C mn between the 
m-th and n-th spins at arbitrary orientation of the magnetic field is performed using the 
software based on the MatLab package which allows us to consider systems of up to ten 
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spins. 



A. Spin system at adiabatic demagnetization 

We consider a linear system of N dip ole- coupling nuclear spins I = 1/2 in an exter- 
nal magnetic field when the Zeeman energy of the order of or even less than the dipolar 
interaction one. The Hamiltonian of the system can be presented in the following form 

H = H Z + H dd (1) 

where the Hamiltonian H z describes the Zeeman interaction between the nuclear spins and 
external magnetic field 

N 

H z = cj J2lk, (2) 



k=l 



U = 7 



He 



is the energy difference between the excited and ground states of an isolated 
spin, Hq is the external magnetic field, 7 is the gyromagnetic ratio of a spin, 7| is the 
projection of the angular spin momentum operator on the z- axes. The Hamiltonian Hdd 
describing dipolar interactions in an external magnetic field 



H « = E (1 - 3 cos 2 9 jk ) lift - ~ (J//,- + IJlt) 



j<k ' " i k 

\ sin 29 jk [e*»>> (J/J+ + J+Jf ) + (i;i k + Jr J*)] - 5 sin 2 [e- 2 ^ fe /+I fe + + e^Hj I~] } 

(3) 

where r^, ^-fc, and are the spherical coordinates of the vector r*^ connecting the j— th 
and /c— th nuclei in a coordinate system with the z-axis along the external magnetic field, H , 
J^and Ij are the raising and lowering spin angular momentum operators of the j-th spin. 
We consider the case when coo ~ D12 = -Q- (here D12 is the dipolar coupling constant for 

r 12 

the nearest spins), and it is necessary to take into account all the terms of the Hamiltonian, 
and not truncate any ones. 

Let us analyze conditions of the adiabatic demagnetization in the spin system. To prevent 
heat exchange between a spin system and a lattice, the evolution time of the considered spin 
system has to be much shorter than the spin-lattice relaxation time, T\. On the other hand, 



decoherence in a spin system refers to how the system loses the quantum coherence features 
and can be characterized by time of the order of T 2 , the lifetime of the free precession signal 
(T 2 << Ti 11, isl) . Decoherence times, T 2 ~ ( r yHi oc )~ 1 ( Hi oc is the local magnetic field 



created by spins) for spin systems at low temperature typically range between nanoseconds 



and seconds (l6j. In contrast to T 2 , T\ takes values in the range from minutes to hours 
[ill , lis) ]. For example, it was found in LiF at T = 2K that T\ = 50 min for 19 F and Ti =15 
hours for Li joj. If to choose the characteristic time of variation of the external field such 
that tch > T2 , the system can be considered as being in thermal equilibrium at every point 
of time and its density matrix is 

P^exJ-m. (4, 



Z ID 



'12 



where f3 is the inverse spin temperature in units of D 12 , (3 = j^., T is the spin temperature 
and Z = Tr |exp ^— j j is the partition function that encodes the statistical properties 
of a system in thermodynamic equilibrium. 



A spin system can be considered in the thermal equilibrium during evolution if {9] 

« iKc- (5) 



dt 



Tr(H 2 ) 

The mean local magnetic field Hi„ r determines by the following expression Hi oc = -\ ™ T d 2 d 

. Our calculations of this field for all 
cases considered below, N = 2 -j- 10, give H\ oc = (0.8 -j- 1.2) D\ 2 , where D12 ~ 10 G [3], 
which leads to the value of the order of few gausses. For example, in LiF calculation gives 
H loc = 7.77 G [9]. With this value and 7 = 2tt x 4005.5 § Eq (5) gives that ^ « 
1.5 x 10 6 f . If we take ^ = 1.5x 10 5 f the field change time from 6000 G to 10 G equels 
t c h — 4 x 10~ 2 s that satisfies the condition of adiabaticity 7\ >> t c h > T 2 ~ 10 -5 s. 

Thermodynamic variables of the system, such as entropy, S and heat capacity C, can be 
expressed in terms of the partition function or its derivatives. For example, the entropy S 
and heat capacity C are given by 

d(pinZ) 

S = kg — (6) 

and 

C = ksfi^, (7) 



respectively. The magnetization M of the spin system can be also expressed in the terms of 
the partition function 



where B is the external magnetic flux density. 

The variation of the parameters at AD is simulated by the following way. At an initial 
condition (3 = (3 in and H = H Qin the entropy is calculated using Eqs. (4) and (6), the 
external magnetic field is decreased and a new temperature is determined from the condition 
S = const. 

Fig. 1 shows the inverse spin temperature versus the external magnetic field. One can 
see that the inverse spin temperature /3 is proportional to ojq at uq » D\i for all the 
spin systems considered here. For uq < D 12 the inverse spin temperature is practically 
independent of ujq (see the inset of Fig.l). The main change of temperature occurs in the 
range of uq/Du from 5 till 1. At ojq/Du > 5 the inverse spin temperature increases 
approximately linearly with a decrease of the magnetic field. In contrast to the dependence 
of the temperature on the external magnetic field, the magnetization is about constant at 
ujq > 5Di 2 and decreases linearly with the magnetic field decrease at uo < D\ 2 (see the inset 
of Fig. 2). At high fields Uq > 3D i2 , the temperature and magnetization became independent 
of the number of spins in the system. 

In the region of high temperatures (low j3), our results coincide with well-known depen- 
dencies obtained in the framework of the high temperature approximation g^- << 1 for 
iV >> 1 (see dark yellow shot dash-dotted lines in Figs. 1 and 2) [j], [lsj. We can see from 



Fig. 2 that the high temperature approximation is valid for uq > -D12 which corresponds to 
< 2.3 (Fig.l). Below we will consider entanglement measures and temperature and field 
dependence of entanglement of the linear chain of dipolar coupled spins under AD. 



B. Generation of entanglement at adiabatic demagnetization 



There are several parameters characterizing the entangled state of a spin system: the von 
Neumann entropy entanglement of formation, log negativity, concurrence of a pair of spins, 



and etc. 



17H20j . We will characterize the entangled states by the concurrence between 



two, m-th and n-th, spins which is defined as 17] 
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C mn = max {q mn , 0} , (9) 

with q mn = \mn - A™ - \mn — \mn- Here A™ = 1,2,3,4) are the square roots of 
eigenvalues, in descending order, of the following non- Hermit ian matrix: 

Rmn = Pmn {? V ® Vy) Pmn {<7y ® CTy) , (10) 

where p mn is the reduced density matrix. For the m-th and n-th spins, the reduced density 
matrix p mn is defined as p mn = Tr mn (p) where Tr mn (...) denotes the trace over the degrees 
of freedom for all spins except the m-th and n-th spins. In Eq. (10) p mn is the complex 

-i 

conjugation of the reduced density matrix p mn and a y is the Pauli matrix a y = . For 

1 

maximally entangled states, the concurrence is C mn = 1 while for separable states C mn = 0. 
The concurrence for two-spin system in the thermal equilibrium is zero at 9 = and 9 = it, 



22 1 . We have confirmed these results 



and it reaches the maximum values at 9 — \ and (p = 
by numerical calculations for all spin pairs in the chains with N = 3, . . . , 10. Therefore, we 
will restrict ourselves to the cases corresponding to the concurrence maximum. 

Fig. 3 presents the results of numerical calculations of the concurrence and heat capacity 
as functions of the magnetic field at AD. With decreasing the magnetic field, the both 
quantities increase up to the maximum and then decrease till zero. Both, the concurrence and 
heat capacity decrease with increasing number of spins for and are practically independent 
of the spin number for . 

Appearance of entanglement at AD and behavior of the concurrence depend on initial 
conditions: initial values of the magnetic field and temperature from which AD is started. 
This dependence is illustrated for a two-spin chain by the phase diagram of Fig. 4. If AD 
starts at point A (Fig. 4) above the curve separating the entangled and separable states, a 
non-zero concurrence is observed till point a\ where the AD phase trajectory crosses the 
boundary line (Figs. 4 and 5). Starting from point B we achieve the entangled state 
between points b\ and 62, while starting from point C does not lead to appearance of the 
entangled state. The similar results are obtained for all neighboring spins in the chains with 
N = 3,..., 10. 

The results for remote spins are presented in Figs. 6-8. At ujq/Di 2 > 5.5, the concurrence 
for the next-neighbor spins C13 has two maximums which are separated by area with zero 
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concurrence (Figs. 6 and 7). The zero concurrence area expands with increasing the spin 
number in a chain (Fig. 6) and the initial temperature (Fig. 7b) and reduces when the initial 
magnetic field increases (Fig. 7a). At AD the first maximum (at higher field) decreases with 
the increase of the spin number in a chain (Fig. 6) and of the initial magnetic field and 
temperature (Fig. 7). The first maximum position moves towards higher magnetic fields with 
increasing spin number and initial temperature. The second maximum also moves towards 
higher magnetic fields with increasing the initial field (Fig. 6) but its position does not 
practically depend on the initial temperature (Fig. 7). 

The concurrence for next-next-neighbor spins C14 also has two maximums while only one 
maximum presents in the magnetic field dependence of the concurrence for the first and fifth 
spins C15 (Fig. 9). In all considered cases, the concurrency for remote spins practically does 
not depend on the number of spins for N > 7. 



C. Discussion and conclusion 



We have investigated entanglement in a linear chain of dipole-coupling spins s = 1/2 at 
AD which is one of the ways to achieve ultra-low spin temperatures. 

The inverse temperature at which entanglement appears is ~ 2.3 . Let us estimate this 
temperature for fluorine in calcium-fluoride CaF 2 with the cipo ar interaction energy of the 
order of a few kHz (in frequency units) Taking as in [l4 , [21 1 Hq = 3 G we have ujq = 12 



kHz, which leads to T = 0.34 //K. The estimated value of temperature is in good agreement 
with those reported early for the spin system s — 1/ 2 with the XY Hamiltonian [4| and 
for the Hamiltonian of two dipolar coupling spins 22]. It is interesting that the transition 



to the ordered states, such as antiferromagnetic, of nuclear spins was observed in a single 
crystal of calcium-fluoride precisely at T = 0.34/iK HQ. 

An unexpected behavior of the next nearest-neighbor and next-next nearest-neighbor 
concurrences was obtained: remote spins in the thermal equilibrium are entangled in two 
distinct temperature and magnetic field regions (Figs 6-8) . In most studied cases t he g round 
zero temperature is entangled [l|, |2j, |23J . However examples were found |20|,|24| where 
the ground state is separable and there are two temperature regions of entangled states. We 
show that the same behavior takes place also for the field dependence of the concurrence 
(Fig. 6) at certain initial values of the temperature and magnetic field (Fig. 4). 



Our investigation demonstrates that the qualitative behavior of entanglement with tem- 
perature and magnetic field can be much more complicated than might otherwise have been 
expected. Thus, the AD technique allows one to generate entangled states between dipolar 
coupling spins in a linear chain. It opens a simple and effective way to the experimental 
testing of entanglement in spin systems. 
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Figure captions 

Fig. 1 Fig. 1 Inverse spins temperature as a function of applied field, , at AD in chains 
with various number of spins: black solid line - N — 4, red dashed line -N = 5, green dotted 
line -N = 6, blue dash-dotted line -N = 7, cyan dash-dot-dotted line -N = 8, magenta 
shot-dashed line - N = 9, yellow shot-dotted - N = 10, dark yellow shot-dash-dotted - the 
high temperature approximation for A^ >> 1. 

Fig. 2 Magnetization, as a function of applied field for an adiabatic demagnetization in 
chain with various number of spins: black solid line - N = 4, red dashed line - N = 5 green 
dotted line -N = 6, blue dash-dotted line N = 7 cyan dash-dot-dotted line -N = 8, 
magenta shot-dashed line - N = 9, yellow shot-dotted A" = 10, dark yellow shot-dash-dotted 

- the high temperature approximation for AT >> 1. 

Fig. 3 Field dependence of concurrence between the first and second spins (a) and of heat 
capacity C per spin (b) in various chains: black solid line -N = 2, red dashed line - N = 3, 
green dotted line - N = 4, blue dot-dashed line - N = 5, cyan dash-dot-dotted line - N = 6, 
magenta shot dashed line - N = 7, yellow shot dotted - N = 8, dark yellow shot dash-dotted 

- JV = 9, navy - JV = 10. 

Fig. 4 The phase diagram for a two-spin system at AD started from various initial points: 
A{ cu = 2.4,(3 = 1.06}; B{ co = 3, (3 = 0.7}; C{u = 2.7, (3 = 0.443}. The black solid line 
is the boundary between the entangled and separable states; green dashed, red dotted, and 
blue dot-dashed lines show the phase trajectories started from A, B, C respectively. 

Fig. 5 Field dependence of the concurrence at AD starting from initial conditions corre- 
sponding to A{co = 2.4, (3 =1.06} (green dashed line) and B{ u = 3, (3 = 0.7} (red dotted 
line) in the Uq — (3 diagram of Fig. 4. 

Fig. 6 Field dependence of concurrence for the first and third spins in various chains at 
AD starting from [3 = 2: black solid line - N = 3, red dashed line -N = 4 green dotted 
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line -N = 5, blue dash-dotted line -N = 6 cyan dash-dot-dotted line -N = 7, magenta shot- 
dashed line - N — 8, yellow shot-dotted N = 9, dark yellow shot-dash-dotted - N — 10. 

Fig. 7 Field dependence of concurrence for the first and third Ci 3 spins in 3-spin chain (a) 
and 6-spin chain (b) at AD starting from various initial points. 

a) - initial inverse temperature (3 = 1 and various initial external magnetic fields: black 
solid line - uq/Di 2 = 5; red dashed line - uiq/Di 2 = 5.3; green dotted line - loq/Di 2 = 5.5; 
blue dash-dotted line - uq/Di 2 = 6. 

b) - initial external magnetic field uj^jDyi = 20 and various initial inverse temperatures: 
black solid line - /3 — 0.3; red dashed line - /3 — 0.4; green dotted line - /3 — 0.6; blue 
dash-dotted line - f3 — 0.9; cyan dash-dot-doted line - f3 — 1; magenta shot-dashed line - (5 
= 3. 

Fig. 8 Field dependence of concurrence Cu for the first and fourth spins in various chains: 
(a) N = 4; (b) black solid line - N — 5, red dashed line - N — 6 green dotted line -N = 7 
for the initial inverse temperature (3 = 2 and magnetic field uiq/D 12 = 40. 

Fig. 9 Field dependence of concurrence C15 for the first and fifth spins in various chains: 
black solid line - N — 5, red dashed line -N = 6, green dotted line -N = 7; for the initial 
inverse temperature (3 = 2 and magnetic field uiq/ D 12 = 40. 
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